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Abstract: In this study, the slope deflection method was presented for structures made of small-scaled 
axially functionally graded beams with a variable cross section within the scope of nonlocal elasticity 
theory. The small-scale effect between individual atoms cannot be neglected when the structures are 
small in size. Therefore, the theory of nonlocal elasticity is used throughout. The stiffness coefficients 
and fixed-end moments are calculated using the method of initial values. With this method, the 
solution of the differential equation system is reduced to the solution of the linear equation system. 
The given transfer matrix is unique and the problem can be easily solved for any end condition and 
loading. In this problem, double integrals occur in terms of the transfer matrix. However, this form is 
not suitable for numerical calculations. With the help of Cauchy’s repeated integration formula, the 
transfer matrix is given in terms of single integrals. The analytical or numerical calculation of single 
integrals is easier than the numerical or analytical calculation of double integrals. It is demonstrated 
that the nonlocal effect plays an important role in the fixed-end moments of small-scaled beams. 


Keywords: slope deflection method; nonlocal elasticity; size effect; axially functionally graded 
materials; transfer matrix; method of initial values 


1. Introduction 


In the present study, static analysis of axially graded nonlocal Euler-Bernoulli beams 
was performed using the slope deflection method. Firstly, the basic equations of a nonlocal 
Euler-Bernoulli beam subjected to distributed load are obtained [1—4]. Then, it is assumed 
that the modulus of elasticity and the moment of inertia functionally change along the axis 
(AFGM) in order to provide a general solution to the problem. For the solution to the prob- 
lem, the transport matrix is formed with the help of Cauchy’s repeated integral. Stiffness 
coefficients were obtained by examining different boundary conditions for different loads. 
Although there are many studies on axially graded beams in previous studies, there is no 
research in the literature on the solution of unit displacement constants of nanoscale frames. 
In the above-mentioned studies, the slope deflection method has been solved for beams 
with constant cross sections and includes a classical structural solution. In the present work, 
the equations are given in the most general form and the solutions can be easily solved 
independently of the changes in the cross section and material parameters. Moreover, as 
this study is examined within the scope of nonlocal elasticity theory for Euler-Bernoulli 
beams, it brings great innovation to the literature in terms of static analysis of both macro- 
and nano-scale systems. 

There are many approaches to solving mechanical problems. Most of these theories 
are methods used in classical elasticity theory. However, this situation differs in the 
analysis of nanoscale structures. It has been repeatedly emphasized that the effect of 
atomic size should be taken into account in the analysis of nanoscale structures. In this 
context, approaches suggesting enriched deformation theories such as nonlocal elasticity 
theory, Reddy’s higher-order deformation theory, and strain gradient theory have been put 
forward by early researchers in analyzing such structures [5]. 


Appl. Sci. 2023, 13, 4814. https: //doi.org/10.3390/app13084814 


https://www.mdpi.com/journal/applsci 


Appl. Sci. 2023, 13, 4814 


2 of 21 


In the literature, there are many methods for solving statically indeterminate systems, 
such as the force method and the slope deflection method. The slope deflection method is 
the most widely used among these approaches, especially for systems with a high degree 
of hyperstaticity. 

Simplifying the application of the slope deflection method, new reduced equations 
based on incorporating both the rotations of the joints and the effects of fixed-end moments 
into one term instead of two terms are presented [6]. The slope deflection method was 
first introduced by George A. Maney to analyze beam and frame systems in the early 20th 
century. This method is basically based on the solution of the continuity equations obtained 
depending on the displacements of the joints in beam and frame systems. The main 
advantage of this method is that it can be easily applied to statically indeterminate systems 
and is suitable for programming [7-9]. Statically indeterminate space frame systems were 
analyzed using the slope deflection method [10]. The slope deflection method was used 
by the researchers to analyze the orthotropic bridge support plates analytically [11]. Static 
analyses of Timoshenko beam—column connections were investigated using the slope 
deflection method [12]. 

With the developing technology, smart materials have replaced the traditional mate- 
rials used in engineering problems. It is of great importance to determine the properties 
of advanced technology materials such as functionally graded materials at a good level in 
order for the structural elements used in almost every field of engineering to have higher 
mechanical strength. 

AFGMs are advanced technology materials used in building elements such as wind 
turbines, helicopter blades, and airframes. Studies have shown that AFGMs subjected to 
loads in different directions show a higher mechanical strength than metal-ceramic-type 
FGMs. However, owing to the fact that the equations of motion of AFGMs can take very 
complex forms, their analytical analysis is quite difficult [13-16]. In the studies conducted 
by the researchers, notably in the last few decades, it has been observed that the resistance 
of axially graded beams to buckling and dynamic behavior is better than that of isotropic 
beams. Researchers analyzed the vibration behavior of AFGM beams using a variation 
iteration approach [17]. The buckling behavior of a composite beam subjected to an axially 
varying load type was investigated using the Ritz method [18]. A group of researchers 
investigated the buckling and vibration behavior of axially graded GPL composite beams 
using the Ritz method. In this context, the mechanical behavior of Timoshenko beams for 
five different load types was examined and, at the end of the study, it was revealed that 
the grading pattern depends not only on the buckling load, but also on the grading type of 
the load (AVLs) [19]. The natural frequencies of a rotating AFGM beam are investigated 
within the scope of nonlocal elasticity theory [20]. At the end of the study, it was seen that 
the grading parameter has a significant effect on the free vibration behavior of the beams. 

Static analysis of an AFGM beam is numerically solved [21]. The mechanical behavior 
of axially graded metal—ceramic microbeams is studied with the help of a higher-order 
theory [22]. The vibration character of an AFGM beam reinforced by CNTs was investigated 
using a differential approximation method. The study revealed that the frequency value 
is very sensitive to the volume fraction of CNTs [23]. The vibration behavior of AFGM 
beams reinforced with GPLs has been investigated within the context of Reddy’s higher- 
order deformation theory [24]. The natural frequency values of AFGM beams reinforced 
with GPLs have been investigated in previous studies [25]. Many structural elements are 
exposed to mechanical effects such as bending, buckling, torsion, vibration, and thermal 
change due to external conditions. In order to minimize such effects, the geometry and 
specific properties of the building elements used should be thoroughly analyzed. 

The buckling behavior of sandwich beams subjected to axially variable load has 
been demonstrated by a series of studies [26]. The dynamic response of a laminated 
composite beam under axially varying loads (AVLs) is studied with the help of a higher- 
order deformation theory [27]. The effect of porosity shape on the buckling and free 
vibration behavior of a functionally graded beam was investigated [28]. The buckling and 


Appl. Sci. 2023, 13, 4814 


3 of 21 


dynamic characteristics of porous beams reinforced with GPLs were investigated, and 
it was reported that porosity shape and axially varying loads have a sensitive effect on 
these behaviors [29]. A new approach is proposed for the free vibration behavior of a 
non-uniform AFGM beam [30]. The modal behavior of a three-dimensional FGM beam is 
obtained using the finite element method [31]. The natural frequencies of the exponential 
FGM beams are presented with an analytical solution method [32]. The free vibration 
behavior of AFGM beams under different boundary conditions has been demonstrated in 
studies [33]. 

Functionally graded materials were first introduced by a group of researchers in 
the Sendai region of Japan and have attracted the interest of researchers [34]. With their 
unique physical and mechanical properties, FGMs have significant potential in the field of 
mechanics in the analysis of building materials. Previous studies have demonstrated the 
mechanical behavior of FGM structures such as bending, buckling, vibration, and thermal 
stress. Investigations have clearly demonstrated the analysis of FGM structures under 
different loading conditions within many studies. In the study presented here, the subject 
was investigated with a different approach related to the work of the researchers and 
briefly discussed [35]. A study conducted by early researchers demonstrated the elasticity 
solution of Euler-Bernoulli beams made of FGM subjected to a transverse load [36]. The 
thermal buckling and vibration behavior of a sandwich plate with a viscoelastic core have 
been studied by taking into account the effect of temperature [37]. Researchers studied 
the thermoelastic behavior of FGM beams by obtaining polynomials of the exact results 
of different equations expressed with interpolation functions [38]. Wave propagation 
behavior of FGM beams under high-frequency loading was analyzed using a finite element 
approach [39]. The ambient-temperature-dependent static behavior of beams in an FGM 
structure was investigated with the help of a high-order beam theory [40]. A new model was 
presented for calculating the buckling load of columns made of FGMs [41]. The mechanical 
behavior of Euler—Bernoulli beams was presented in which material properties vary in 
either the horizontal or vertical direction [42]. Functionally graded materials are a type of 
composite material whose material properties are designed to show continuity from one 
surface to another. As for laminated structures, stress discontinuities occur between two 
surfaces. These discontinuities are eliminated in FGMs. Owing to these properties, FGMs 
are widely used today in many areas such as the design of medical devices, as well as in the 
fields of construction and ship and aircraft engineering. Beams made of FGMs are classified 
into two basic groups according to their design status and their grading according to the 
thickness and length of the element, with their mechanical properties varying in different 
directions (axially or longitudinal) [43,44]. These types of beams are called axially graded 
beams (AFGs). There are a limited number of studies in the literature on the analytical 
approach of AFGs [45,46]. There are many methods in the literature to determine the natural 
frequency of non-uniform axial FGM beams. The bending and buckling behaviors of bi- 
directional FGM beams are investigated within the scope of strain gradient theory [47,48]. 
These are the complementary functions approach [49], finite element method [42-50], 
power series approach [51] used to define mode shapes, Chebyshev polynomials [52], 
differential quadrature method [53], differential form including initial values [54], and 
Fredholm integral equations [55]. Classical elasticity theory is not sufficient to analyze 
the behavior of nanoscale structures in the most accurate way. Therefore, it has emerged 
that an enriched continuum model is needed to describe the mechanical behavior of such 
structures. Although the molecular dynamics approach (MD) is a method used to analyze 
such structures, it contains a high level of computational load, making it difficult for such 
simulations to express a general solution. The nonlocal elasticity theory approach states that 
the stress field at a point is affected not only by the strain at that point, but also by the strain 
of all adjacent points. In this context, nonlocal elasticity theory has an important place in the 
analysis of size-dependent mechanical behavior [56,57]. In recent years, numerous studies 
are presented to display the importance of the size effect in small-scale structures [58-60]. 
The slope deflection method is a method based on expressing the displacement components 
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at the nodes as unknowns. There are many publications in the literature related to this 
method. However, there is no publication describing the method in nonlocal theory for 
structures composed of AFG elements. Stiffness coefficients and end moments are analyzed 
within the scope of nonlocal elasticity theory using the method of initial values (MIV). 


2. Basic Equations 
The basic Equation (1) [1,2,5] of a nonlocal Euler—Bernoulli beam is as 


dv = dp _ = =M-’M" dM _, dT _ 
dz © dz E(z)I(z) ’ dz " dz p(2) v) 


Note that v, gy, M, T, E(z), I(z), and ¥ are deflection, rotation, moment, shear force, 
modulus of elasticity, moment of inertia, and small-scale parameter revealing a small-scale 
effect on the response of structures in nano size, respectively. Here, p(z) is the distributed 
load [61-68]. Using Equation (1), the above expression can be written as 


6) 0 1 0 0 6) 
dl o|_|9 0 -zary 1 @ — pe 
az | M 0 0 , i||™m + : (z) (2) 
T 0 0 0 0 T —p(z) 
Equation (2) can be expressed as 
dy ,.-,2 > 
dz =4Y a f; ( ) 0 (3) 
The solution of the MIV given by Equation (3) is 
> => Zz = 
y(z) = TM(2).¥o+ f TM(z—T).f (t)dt (4) 


where T is parameter used on variables in double integrals, « is a variable that shows the 
variation of E and I along the axis, and TM(z) is the transfer matrix. The transfer matrix is 


obtained for variable E and I using Cauchy’s repeated integration formula as below [63]. 


Lz So wijity’® —Jo Fees 

TM(z)= |° 1 —Jo Err’ Jo raya (5) 
0 0 1 Zz 
0 0 0 1 


3. Applications of the Initial Values Method 
3.1. A Beam Subjected to a Distributed Load 


Mo, To, Yo, and po represent the initial values of the bending moment, shear force, 
deflection, and rotation, respectively. The unknowns of a beam subjected to a distributed 
load are calculated using Equations (4) and (5) as below (See Figure 1a). 


0(2) = —Mo(Jo asia) + Jo PO (Io "Staite + Ree Jat 
To(Jo ratey4t) + 0+ goz 
p(z) = —Mo( Jo ranay@*) — To (Jo ratnay@*) + 99 
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Figure 1. (a) A beam subjected to a distributed load; (b) a beam subjected to two different loads. 


3.2. Beam Subjected to Two Different Distributed Loads 
In this case, Equation (4) can be written for both regions as follows (see Figure 1b). 


Vi (2) = TMG@).Ho + [TM -7).f (r)dz, 0 <2 < qh (9) 
Vo(z) = TMG@).H,(9L) + [TM —1).fa(r)dz,0<Sz< (1-0) 
where 
0 0 
=s _ pil) = _ ¥p2(z) 
fi) =| FENG) |, f,(2) = | Fane) (11) 
—p1(z) —p2(z) 


The unknowns of a beam subjected to two distributed loads are calculated using 
Equations (9) and (10) as below (see Figure 1b). 
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3.3. Beam Subjected to Singular Force 
In this case, p(z) = 0, Equation (4) can be written for both parts as follows (see 


Figure 2a). 
¥1(2) =TM(z).Vo, 0 <2 < HL (19) 
Voz) = TM(z).(y4(7L) + B), 0<2< (1-W)L (20) 
anaes 
= (0 0 0 -@)" (21) 
Q p(zZ) 


si nLE@ UZ) (-n)L E(z) (z) 
——3#764»7<+—- 
er 
(a) (b) 


Figure 2. (a) Beam subjected to singular force; (b) both ends of a fixed beam subjected to a distributed load. 


The unknowns of the problem can be written as follows. 
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4. Calculation of Fixed-End Moments 
4.1. Fixed Beams with Both Ends Subjected to the Distributed Loads 


The positive directions for the end moments are different in MIV and SDM (see Figure 3). 


M. M.. M.., 
G4 $y ¢ 2 
i j MM J 
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Figure 3. (a) Positive directions in MIV; (b) positive directions in SDM. 


Appl. Sci. 2023, 13, 4814 7 of 21 


A fixed beam with both ends subjected to the distributed load is shown in Figure 2b. 
In this case, two of the initial values and end conditions can be written as 


v9 = 0, Yo =0; v(L) =0, g(L) =0 (28) 


The other two initial values are found with the help of Equations (6) and (7) as follows. 


Mo = ((Jo° grainy) * Jo PO (So * “Stamey + Foon) 47 
— (Jo etary 4#) * Jo PCO (Jo Faye rare) 47) 2 
i arty) Jo Eratray 4 — (Jo zrajnayte) Jo’ eeasray@) 
a resting ®) Jo PC) (Jo * exafnay — aera) at 
“( sayrayt) * Jo P(t) (So —seartay ae + Beret) 4) (30) 
/ (i sts ) Jo wera — (Jo earey@®) Jo asters) 


Special cases 


(a) Uniform beam with constant cross section subjected to the uniform load (see Figure 4a). 
In this case, 


p(z) =p, Ez) =E, Itz) =1 (31) 
Using Equations (29) and (30), the initial values Mg and To are found as follows. 
1 
Mo = —=<p (127 +L’), To = (32) 
The variation of the initial moment (Mpo) in classical and nonlocal theory with respect 
to the 7/L ratio is given in Figure 5a. 


(b) Uniform beam with a constant cross section subjected to the triangular load (see 
Figure 4b). In this case, 


p(z) = £2, E(@) =E, 12) =1 (33) 
Using Equations (29) and (30), the initial values Mp and To are found as follows. 


2 
L*p Beep a (34) 
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The end moments are found using Equation (8) as below. 


oe TT 
M(L) = —-p(207" +17), T(L) = -T? - —? (35) 


The variation of the end moment (M_) in classical and nonlocal theory with respect to 
the } ratio is given in Figure 5b. 
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Figure 4. (a) Beam subjected to the uniform load; (b) beam subjected to the triangular load. 
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Figure 5. (a) The variation of the initial moment (Mo); (b) the variation of the end moment (Mr). 
4.2. Fixed Beams with Both Ends Subjected to the Different Distributed Loads 


A fixed beam with both ends subjected to the two distributed loads is shown in 
Figure 6. In this case, the initial values and end conditions are 


V9 = 0, po =0; v2((1—)L) =0, g2((1—4)L) =0 (36) 


Pw) P(2) 


| (I-n)L 
mm NL (1-n) 


—_——_ 
Figure 6. A fixed beam with both ends subjected to the two different distributed loads. 


The other two initial values are found with the help of Equation (15) as follows. 
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Special cases: A uniform beam subjected to different types of loads. 


(a) Uniform beam with a constant cross section subjected to two different uniform loads 
(see Figure 7a). 


For 7 = 5 and 
PiZ)=pix USe< HL 


39 
p2(z)=po, O<z<(1—y)L Ma 
Mo and M2(L/2) are obtained using Equations (17) and (37) as follows. 
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The variation of the initial moment (Mp) in classical and nonlocal theory with respect 
to the 7/L ratio is given in Figure 8a. 
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(b) Uniform beam with a constant cross section subjected to the two different triangular 


loads (see Figure 7b). 
For 7 = 7 
pi(z) = "BZ, 0<z<yL 
po(z) = pr—- =, O<z<(1—7)L (42) 
Mo and M>(L/2) are obtained using Equations (17) and (37) as follows. 
Mo = Jp (—16p1 (157° + L*) — 9L* po) (43) 
1 
Mp(L/2) = gen (—16 po(157° + L*) — 9L* pi) (44) 
Tp and T>(L/2) are found using Equations (18) and (38) as below. 
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The variation of the initial moment (Mpo) in classical and nonlocal theory with respect 
to the 7/L ratio is given in Figure 8b. 
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Figure 7. (a) Uniform beam with a constant cross section subjected to the two different uniform loads; 
(b) uniform beam with a constant cross section subjected to the two different triangular loads. 
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Figure 8. (a) The variation of the initial moment (Mo); (b) the variation of the initial moment (Mo). 


4.3. Beams Subjected to Singular Force 
4.3.1. Fixed Beam with Both Ends Subjected to Singular Forces 


In this case (see Figure 9a), two of the initial values and end conditions are 


%=0, go =0; va((l1—4)L)=0, — ga((1—4)L) =0 (46) 
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Figure 9. (a) Fixed beam with both ends subjected to singular force; (b) fixed-pinned beams subjected 
to singular force. 


Two other initial values (Mg and Ty) are obtained using Equations (25) and (26) 


as follows. 
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Special Case: E and I are constant. 
For E(z) = Eand I(z) = I, Mo and Tp are found with the help of Equations (47) and (48) 
as follows. 
Mo = (q—1)?9(-L)Q, To = (q—1)?(2y +1) Q (49) 
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4.3.2. Fixed-Pinned Beams 


In this case, two of the initial values are (see Figure 9b) 
%=90, go=9 (50) 
The end conditions are 
v2((1—y)L)=0,  Mo((1—y)L) =0 (51) 


Mp and To are obtained using Equations (25) and (26) as below. 
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Special Case: E and I are constant. 
For E(z) = Eand I(z) = I, Mo and Ty are found with the help of Equations (52) and (53) 
as follows. 


My =—5 (1-2) (q-1)9LQ, Ty = 5((4-3) 1? +2) Q (64) 


5. Calculation of Stiffness Factors 
5.1. Calculating the Moment jt Corresponding to the Unit Rotation for the Figure 10a 


In this case, the applied moment p rotates the end j by an amount ¢(L). Two of the 
initial values are (see Figure 10a) 


09 = 0, fo = 0 (55) 


—— 
(b) 


Figure 10. (a) The applied moment p/ rotates the end j by an amount of ¢(L); (b) the applied moment 
y rotates the end j by an amount of p(L). 
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Using Equation (4), the unknowns can be written as 


(2) = Mo ( a Fala Tw) (fa Ea)a) + da) Pe) 


z Zz 

2) = Mo(— f° ararray 4) -(f arayreay a Oa 
M(z) = Mo + Toz (58) 
T(z) =T (59) 

The end conditions are 
v(L) =0, M(L) = hu (60) 
The other initial values (Mo, Tp) are obtained using Equations (56) and (58) as follows. 
My - L —s nate da” We _ Ein eL) 


Sy easy Jo equity’ 
The moment p/ is found with the help of g(L) = —1 as below. 
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Special case: E and J are constant. 
For B(2)-— &,. 1a) =I, 
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12’ 


(63) 


5.2. Calculating the Moment jt Corresponding to the Unit Rotation for the Figure 10b 


In this case, the applied moment p/ rotates the end j by an amount ¢(L). The initial 
values and end conditions are (see Figure 10b) 


v9 = 0, My =0; o(L) =0, M(L) =p (64) 


Using Equation (4), the unknowns can be written as 
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M(z) = Toz, T(z) = To (66) 


The other initial values (@o, Tp) are obtained using Equation (66) as follows. 


L a«(L—«) 
_v Jo E(aylay 4 T = 67) 
D ’ L 
The moment pi and Tp are found with the help of g(L) = —1 as below. 
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Special case: E and I are constant. 
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5.3. Calculating the Moment jt Corresponding to the Unit Deflection for the Figure 11a 


In this case, vertical displacement occurs at the j-end of the beam owing to the moment 
yu. The initial values and end conditions are (see Figure 11a) 


v9 = 0, po =0; g(L) =0, M(L) = p (71) 


E(z) (2) 
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Figure 11. (a) Vertical displacement occurs at the j-end of the beam owing to the moment p; (b) 


vertical displacement occurs at the j-end of the beam owing to the moment pi. 


Using Equation (4), the unknowns can be written as 


v(z) = Mo ( ‘ Faia) in) To ( a in) (72) 
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M(z) = Mo + Toz, T(z) = To (74) 


The other initial values (Mo, Ty) are obtained using Equations (73) and (74) as follows. 
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The moment pi is found with the help of v(L)= 1 as below. 
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5.4. Calculating the Moment y Corresponding to the Unit Deflection for the Figure 11b 


In this case, vertical displacement occurs at the j-end of the beam owing to the moment. 
The initial values and end conditions are (see Figure 11b) 


v9 =0, Mo = 0, M(L) = pn, g(L) =0 (78) 


Using Equation (4), the unknowns can be written as 
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M(z) = Toz, T(z) = To (80) 
The other initial values (@, To) are obtained using Equation (80) as follows. 
La 
B So Earay 
99 = =F (61) 
The moment pi is found with the help of v(L) = 1 as below. 
L 
ne L L a(L—«a) (82) 
L fo Emir t® — Jo Fayre 44 
Special case: E and I are constant. 
Por e(gy=£,1(2)—1, 
3EI 3EI 
To = T(L) = Fy (83) 


neo [2 ’ 


6. Examples 
6.1. A Frame Subjected to Triangular Load 
In this example, the moments at each joint of the frame shown in Figure 12a are 


obtained in the frame of nonlocal elasticity. 
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Figure 12. (a) The frame subjected to the two triangular loads (p); (b) the variation of rotation (6g) at 


point B. 
Equilibrium equations can be written in nodes B and C as follows [64] (see 


Equations (43) and (63)). 
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ron On + (SEE FP) 6c — Sp (4877 +5127) =0 (85) 
Solving Equations (84) and (85) simultaneously, 0g and 6c are found as follows. 
Ly Lop (4877 + 5L2”) nea (487711 Lz + 5L1L23) (86) 
oe 192EI (B Ly + 2L2) 


o3 = — DET (BL + 2Ly) 


Appl. Sci. 2023, 13, 4814 16 of 21 


In the limit y goes to zero, 6g and @¢ are 


lim 63 = SL Lop lim 6¢ = 5L1 Leip 
con 192BEIL, +384EIL2’ 70 © 192BEIL, + 384EIL2 


(87) 


The unit displacements obtained as a result of the present work are validated by the 
results given in [64] for L] = 12m, L7 = 8 m, B = 1, y = 0, and p = 24 kN/m, for 03 and 6¢ 


(see Table 1). 
960 960 


0 — 
Tel’? “FEI 
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Table 1. Unit displacements for a frame system subjected to vertical triangular load [64]. 


Unit Displacements Reference Study (7/L = 0) [64] Present (7/L = 0) 
960 960 
OB FFT ~~ TEI 
960 960 
9c FET TEI 


Special Case: Variable cross section. 
In the frame shown in Figure 12a, the variation of the moment of inertia and modulus 
of elasticity function of the elements AB and DC are (see Figure 13) 


Ce ene ees ee ee (89) 
Ly Ly 
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Figure 13. The variation of the moment of inertia of the elements AB and DC (see Figure 12a). 
For L; = 12 m, mj 9, is calculated using Equation (62) as follows. 


Equations (84) and (85) take the following form: 
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6g and @¢ are obtained by solving the system given in Equation (91). 
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8 EI (33.1673B2 + 13.3001BL2 + Lz?) . 
For Lj = 12m, L2 =8m, B=1, 7 =0, and p = 24kN/m, 
_ 93.9518 (03) 
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6.2. A Frame Subjected to Horizontal Triangular Load 


In this section, a frame system subjected to non-symmetrical load is given as shown in 
Figure 14. Equilibrium equations can be written in nodes D and C and as follows [64] (see 
Equations (34), (35), (63) and (77)). 


4EI , 4BEL | 2BEI | 2BEIOc , 6EIS — PLY _ 
(get eta tere, 4 a a (94) 
4EI , 4PEI | 2BEI0p , GEIS _ 
(4 Te) Oc 4 ora 3 0 
| | 
| El, Li EI, Li | 
| | 
P A B 
Figure 14. Frame subjected to non-symmetrical loading (p). 
The equilibrium condition in the AB direction is 
4EI 2EI 4EI 2EI 
eo iy Ly Ly 2° 20° 


6c, 8p, and 6 are found by solving Equations (94) and (95) as follows. 


ply Lz (60By7Ly —17BL} +12077L2—12L3 Ly ) 


Ac 240EI(6BL1+L2)(BL1+2L2) 
pLy Lz (60Py7Ly +7BL3+12077L7—8L4L> ) 
0p = (96) 
240E1(6BL1+L2)(BLy+2L2) 
j= pL? (—606-77L1+9BL3—4077L2 +417 Ly ) 


240ET(6BL4 +L) 


The unit displacements obtained as a result of the present work are validated by the 
results given in [64] for L} =3m, Lp = 3m, B = 1, y =0,and p = 30 kN/m, for 6c, 8p, and 6 
(see Table 2). 


Table 2. Unit displacements for a frame system subjected to horizontal triangular load [64]. 


Unit Displacements Reference Study (7/L = 0) [64] Present (7/L = 0) 
4.66 4.66 
Ac ET ET 
0.161 0.161 
OD ~ EI ~ EI 
6 18.8 18.8 
T EI 


The end moments can be written as below. 
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2EI 4EI 
_ | Q 
Mpc ( iz 9c +4 iz D) (99) 
2EI 4EI 
Neg=- Oo be (100) 
Li i 
2EI 6EI 
== aay 101 
Mec i 6c+ B (101) 
4EI 6EI 
— | 1 2 
Mcp (¢ 9c + re) (102) 


In the limit 7 goes to zero, the end moments (see Table 3) revert to the moments given 
in [64]. 
Table 3. End moments for a frame system subjected to horizontal triangular load [64]. 

End Moments 


Reference Study (7/L = 0) (kN.m) [64] Present (7/L = 0) (kN.m) 


Mis 25.9 —25.9 
Mpa 3.32 3.32 
Mpc 3.39 9.52 
Men —6.32 =655 
Mac 9.43 9.48 
Mcp —6.32 =637 


The variation of the rotation at the joint D in classical and nonlocal theory with respect 
to the 7/L ratio is given in Figure 15. 


El@p 


Bp 


a 
0.06 0.07 0.08 0.09 0.10 £2 
Figure 15. The variation of the rotation at the joint D. 


7. Conclusions 


In this paper, the slope deflection method for AFG beams with a variable cross section 
was presented in the frame of nonlocal elasticity. First, the transfer matrix was presented for 
AFG beams with a variable cross section. The stiffness coefficients and fixed-end moments 
were calculated systematically using MIV. Fixed-end moments of a beam subjected to 
distributed loads include the nonlocal parameter (7) (see Section 4). It was observed that 
the nonlocal parameter was not included in the stiffness coefficients. That is, stiffness 
coefficients are independent of the nonlocality parameter (see Section 5). The difference 
between the moment values obtained by local and nonlocal theories grows significantly 
when the 7/L ratio increases (see Figures 5a,b and 8a,b). A similar situation also occurred 
in the rotation of joints in the frames (see Figure 12b). Therefore, the fixed-end moments 
calculated by nonlocal elasticity should be used in the static analysis of the nanostructures. 
All of the obtained results were validated by comparison to the reference study. 
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With the help of the present work, the unit displacement constants of the axially 
graded beams in different loading conditions were calculated using the given formulas. 
The obtained results were validated by comparison to the reference study. As statically 
indeterminate systems may contain many unknowns, solving such systems statically 
increases the computational load. This study offers an innovative approach in terms 
of easier analysis of systems such as frames of different dimensions that we encounter 
in engineering structures. The fact that the 7/L ratio is zero in the study indicates the 
abandonment of the nanoscale and the transition to the classical state. Another important 
aspect of this study is that the relevant solutions are in a general form. That is, the solutions 
do not depend on the dimensions of the problem. 

The authors agree that this study will be an important resource for structural engi- 
neering problems, wind turbine design, relatively small-sized electronic devices such as 
NEMS and MEMS, and the design of biomedical devices. We also believe that a more ad- 
vanced composite material model can be obtained by designing the beam as a bi-directional 
structure, in addition to the structure of axially graded beams in the current study. 
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